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ABSTRACT 

Developing accurate analysis techniques to combine various probes of cosmology is essen- 
tial to tighten constraints on cosmological parameters and to check for inconsistencies in our 
model of the Universe. 

In this paper we develop a joint analysis framework for six different second-order statistics 
of tracers of the dark matter density field, namely galaxy position, shear, and magnification. 
We extend a data compression scheme developed in the context of shear-shear statistics (the 
so-called COSEBIs) to the other five second-order statistics, thereby significantly reducing 
the number of data points in the joint data vector 

We simulate a likelihood analysis for the Dark Energy Survey (DES) in a five dimensional 
cosmological parameter space (plus basic parametrization for bias and correlation parameter) 
comparing the information content of the individual probes to several combined probes (CP) 
data vectors. Given the significant correlations of these second-order statistics we model all 
cross terms in the covariance matrix; furthermore we go beyond the Gaussian covariance ap- 
proximation and use the halo model to include higher order correlations of the density field. 
We find that modeling both, Non-Gaussianity and cross terms, is essential for accurate like- 
lihood contours from the CP data vector; assuming Gaussian errors or uncorrelated probes 
impacts the contours at a similar level as assuming perfect knowledge of galaxy bias. 
We also identify several nuUtests based on the degeneracy of magnification and shear statis- 
tics which can be used to quantify the contamination of data sets by astrophysical systematics 
and/or caUbration issues. 
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1 INTRODUCTION 

High quality data sets from near-term wide-field imaging surveys, 
e.g. Kilo-Degree Survey (KiDS'), Hyper Suprime Cam (HSC-), 
Dark Energy Survey (DES''), allow for tight constraints on cosmo- 
logical parameters from the Large Scale Structure (LSS) of the Uni- 
verse, being complementary with Cosmic Microwave Background 
(CMB) constraints from the Wilkinson Microwave Anisotropy 
Probe (WMAP'') (see Hinshaw et al. 2012, and references therein) 
and Planck^. The improved data quality and the small statistical 
uncertainties (as a result of the increased survey volume) pose new 
challenges for the data analysis; the development and refinement 
of LSS data analysis methods is crucial for the success of even 
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larger, future data sets from the Large Synoptic Survey Telescope 
(LSST*), and from future satellite missions Euclid' and the Wide- 
Field Infrared Survey Telescope (WFIRST**). 
For all of the aforementioned surveys the tightest constraints on 
cosmology will be obtained from a joint analysis of all probes 
that can be extracted from the data (e.g., cluster mass function, 
shear peak statistics, BAO peak fitting, and various second-order 
statistics derived from clustering, shear, and magnification). Com- 
bining LSS with Supemovae and CMB constraints is straightfor- 
ward; due to the fact that these probes have very little correlation a 
joint likelihood analysis frequently comes down to multiplying the 
con'esponding posteriors probabilities (e.g., Kilbinger et al. 2012). 
Combining the various probes of LSS themselves is complicated 
for several reasons: First, the cosmological information of various 
LSS probes is highly coixelated, which prohibits a joint analysis 
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on the level of posterior probabilities. Instead the analysis requires 
a joint likelihood using a covariance matrix that includes all cross 
correlation terms between the individual probes. Second, not only 
is the cosmological information coiTelated, even more problematic 
are the correlations of various systematic effects originating from 
astrophysics and the measurements themselves. 
As a consequence we do not only have to include the covariance 
matrix's cross terms in the likelihood analysis, moreover a joint 
analysis requires a joint modeling of the CP data vector, or more 
precisely, it requires a prediction code for the combined data vector 
that ensures consistent modeling of the considered quantities as a 
function of cosmology and also as a function of the uncertainties 
in the nuisance parameters. Developing such a CP prediction code 
is challenging given that modeling the individual probes already 
requires refined knowledge and high-level expertise on the corre- 
sponding astrophysics and systematics. Although this knowledge 
is present in the corresponding communities, even the individual 
analysis methods are under constant development in order to meet 
the new data quality, and unfortunately these methods are largely 
independent from each other. Phrasing the problem differently: the 
large correlation of the LSS probes is not reflected in the correla- 
tion of the development of the individual analysis techniques. 
For example, probably the most important astrophysical uncer- 
tainty for clustering based measurements is the relation of dark and 
luminous matter, modeled through various bias parametrizations 
and/or Halo Occupation Density (HOD) models. Constraints on 
these models come from measuring cross correlations of shear and 
clustering (sometimes called galaxy-galaxy lensing). Cosmic shear 
uses the same cross terms to offset uncertainties due to intrinsic 
alignment; simply combining both methods uses the galaxy-galaxy 
lensing information twice. Similar problems occur when modeling 
shear calibration which affects cluster masses calibrated through 
weak lensing, cosmic shear, galaxy-galaxy lensing, and shear peak 
statistics, all at the same time. Other examples are the modeling 
of baryonic uncertainties and photo-z calibration which affect all 
probes but in different ways. 

Whereas solving all these problems is beyond the scope of this pa- 
per, it is our intention nevertheless to take first steps towards a co- 
herent analysis framework of LSS probes. We limit our problem 
to second-order statistics (power spectra, correlation functions or 
linear transformations thereof) that can be derived from a galaxy 
catalog containing measurements of galaxy shear, galaxy position 
and galaxy magnification. As we explain further in Sect. 2 we ob- 
tain six different second-order statistics from these measurements 
corresponding to six different probes of the density field. We ex- 
clude galaxy clusters and shear peak statistics for now, since these 
are first-order number count statistics which cause additional com- 
plications in the sense that they require a different likelihood func- 
tion (Poisson distribution instead of Gaussian) in their analysis. In- 
corporating these probes at the level of the covariance matrix in a 
Gaussian likelihood together with second-order statistics is ques- 
tionable. 

In terms of nuisance parameters we restrict our analysis to a simple 
two parameter bias model for probes that involve clustering (see 
Sect. 5.1), however the extension to other astrophysical contami- 
nations, e.g. intrinsic alignment and baiyonic effects is straightfor- 
ward, and although not being part of the analysis we address it in 
the discussion. 

Within the aforementioned restrictions our analysis advances the 
existing state of the art of simulated likelihood analysis: 1) We sim- 
ulate an actual likelihood analysis in a five dimensional cosmolog- 
ical parameter space (and two nuisance parameters), 2) we use a 



full Non-Gaussian covariance which includes correlations between 
all probes and account for higher-order correlations of the density 
field (see Sect. 4), and 3) we develop a data compression scheme 
for the joint likelihood analysis which simultaneously solves the 
shear-shear EB-mode problem. 

This data compression scheme was developed by Schneider, Eifler 
& Krause (2010) (hereafter SEKIO) to solve the problem of cal- 
culating the a shear E-mode two-point statistics (which contains 
the cosmological information) from a given shear-shear correlation 
function on a finite interval. Since information from shear data is 
limited to angular scales [)?miii; J^max] any E/B-mode statistic which 
requires information on larger or smaller scales suffers from so- 
called E/B-mode mixing or leakage. This problem is examined in 
Kilbinger, Schneider & Eifler (2006) for configuration space quan- 
tities, finding a significant (scale-dependent) bias for formerly used 
shear statistics, e.g. aperture mass dispersion or shear dispersion. 
The issue has been addressed in even greater detail for Fourier 
space quantities, mostly in the context of CMB polarization experi- 
ments; several groups developed and refined a Pseudo-Cl technique 
(Hivon et al. 2002; Brown, Castro & Taylor 2005) that has been ap- 
plied to simulated shear data in Hikage et al. (2011). In Fourier 
space E/B-leakage largely depends on the mask of the survey; sev- 
eral mitigation schemes have been developed (Lewis 2003; Smith 
2006; Kim & Naselsky 2010). 

Except for shear-shear none of the other five second-order statistics 
suffers from the EB-mode problem; nevertheless the data compres- 
sion aspects of the COSEBIs are highly desirable for these probes 
as well. Furthermore, the extension of the COSEBIs scheme allows 
for a joint cosmological analysis that involves a clean separation of 
the cosmic shear signal into E- and B-modes. 



2 BASIC CONCEPTS 

We consider the observables shear y, magnification ji, and 
galaxy position g. From these observables the following second- 
order statistics can be obtained: shear-shear (77), magnification- 
magnification (j-ij-i), galaxy position-galaxy position (gg), shear- 
magnification (7^), shear-position {yg), and magnification-position 
(yug). We want to comprise this second-order cosmological informa- 
tion into a COSEBIs data vector 



E = (E'''', W\ E««, E"", E''", E"'')' , 



(1) 



where each E" contains five COSEBI modes (see SEKIO, Eifler 
2011; Asgari, Schneider & Simon 2012, for justification of the 
number of modes). The goal of this paper is to simulate a multidi- 
mensional likelihood analysis, where "simulated" means that E is 
computed from a fiducial cosmological model (see Table 1) using 
our prediction code; we will refer to this data vector as the fiducial 
data vector from now on. 

We assume that the errors of the input data vector E are described 
by a multivariate Gaussian 



i(E|p„) : 



1 



(27r)~/2 VlCi 



exp 



1 



• M)' C"' (E - M) 



(2) 



where p^o denotes the cosmological parameter vector that is as- 
sumed in the model vector M, hence M = M(pco). 
The posterior probability in cosmological parameter space is ob- 
tained via Bayes' theorem 



P(Pco)L(E|p,„) 
P(E) 



(3) 
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where P(Pco) denotes the prior probabiUty (we assume non- 
informative priors) and the evidence P(E) can be calculated as an 
integral over the likelihood P(E) = Jdpco/'(Pco) ^'(E|pco) provid- 
ing a normalization constant for the posterior probability. 
Given the functional form of the likelihood as in Eq. (2) the error 
bars are fully determined by the covariance of the COSEBIs' data 
vector, which correspondingly to the definition in Eq. (1) reads 
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(4) 



with C being symmetric. While postponing a detailed description 
of the covariance's modeling to Sect. 4, we note that we assume 
the covariance to be constant with respect to the point in parameter 
space where the likelihood is evaluated. As shown in Eifler, Schnei- 
der & Hartlap (2009) (for Gaussian shear-shear covariances) this 
assumption is problematic and, depending on the survey parame- 
ters, can have significant impact on the parameter constraints. We 
acknowledge that the covariance matrix, since predicted from a cos- 
mological model, in principle has to vary with respect to cosmology 
(see fee et al. 2012; Kilbinger et al. 2012, for corresponding appli- 
cation to shear data) and we will pursue a corresponding extension 
of this work in the future. 

In practice the COSEBIs are calculated from the correlation func- 
tions of the three observables. The corresponding power spectra are 
related to these correlation functions as 



^""(l?) 
^««(l?) 
^''"(l?) 
^«(l?) 
^''«(1?) 



^ j defjomc^-ie), 
— f dee j2(M)c^''ie), 

In J 

J defjiimc^e). 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



where we point out the J2 in the polar-scalar correlation functions 
and (see e.g. Bartelmann & Schneider 2001, for a deriva- 
tion). We will return to the filter functions in Sect. 5.2. 



3 MODELING THE DATA VECTOR 

In this section we describe the prediction code for the model 
data vector (please see Fig. 1 for an illustration), which is an ex- 
tended version of the shear-shear prediction code described in Ei- 
fler (2011). All projected quantities are computed from the non- 
linear density power spectrum which we calculate from an initial 
power spectrum Pss(k) oc k^^ using the transfer function of Eisen- 
stein & Hu (1998). In order to model the non-linear evolution of the 
density field we develop a Hybrid approach combining information 
from Halofit (Smith et al. 2003) and the Coyote Universe Emulator 
(Lawrence et al. 2010). The latter emulates P^^ over the range k e 
[0.002; 3.4] h/Mpc within z e [0; 1] to an accuracy of 1% for cos- 
mologies within fi,^!^ g [0.120;0.155], Qb/i^ e [0.015; 0.0235], 
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Figure 1. Schematic illustration of the modeling of the CP COSEBIs data 
vector for a given cosmology. 



n, e [0.85; 1.05], erg e [0.6; 0.9], wq e [-1.3; -0.7]. For any cos- 
mology, z, and k within the aforementioned range, we solely rely 
on the output of the emulator. For all other parameters we compute 
the non-linear part of Pss from Halofit and rescale this solution by 
a factor 



/(fc,z,p,.o) 



[jCoyote 



"(fc,z,P^!,'"') 



(11) 



where p^o is the cosmology parameter vector of interest and pf°'"' 
is the closest point in parameter space where the Emulator returns a 
solution (close is defined as minimum difference in each parameter 
separately). 

In order to simulate wCDM models we follow the strategy out- 
lined in icosMO (Refregier et al. 2011), which interpolates Halofit 
between flat and open cosmological models to mimic Quintessence 
cosmologies (please also see Schrabback et al. 2010, for more de- 
tails). Outside the parameter range of the Emulator the precision of 
Pss will of course be significantly below 1%; we nevertheless be- 
lieve that our approach supersedes other implementations of mod- 
eling non-linear structure growth for multiple cosmologies. For ex- 
ample, when using Halofit alone it has been shown that the lensing 
power spectrum is substantially underestimated (e.g., Hilbert et al. 
2009). Throughout this paper we assume a redshift distribution as 
expected fro DES. More precisely, this is modeled by modifying 
a redshift distribution measured from the Canada-France-Hawaii 
Telescope Legacy Survey (see Benjamin et al. 2007, adjusted for 
the lower mean redshift of DES). The exact parameterization reads 



n(z) = A' — exp 



(12) 



with a = 2.0,/? = 1.0, zo = 0.5. 



3.1 Modeling the projected power spectra 

Kaiser (1992, 1998) show that projected power spectra are related 
to the 3D power spectrum of density fluctuations P^s via a Fourier 
equivalent of Limber's equation, i.e. 



Pn(e) : 



P66 (k,x) 



(13) 
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with qi,q2 being weight functions, k = tlx, and fk(x) being the 
comoving angular diameter distance which corresponds to the co- 
moving coordinate x for the case of vanishing curvature. For sim- 
plicity, we will assume the latter in our analysis; note that the tools 
described in this paper are nevertheless independent of this assump- 
tion. 

In case of the shear the weight functions q read 



- f'dx'p.ix')^^, (14) 



where a{x) is the scale factor and pj(x)dx = p,(z)dz is the redshift 
distribution of source galaxies in the ;* tomography bin. We do 
not consider tomography in this paper and drop the corresponding 
denotation of different redshift bins from now on. 
Using these weight functions the expression for the shear power 
spectrum reads 



(15) 



In the weak lensing approximation the magnification yu equals twice 
the convergence k, where the latter equals the shear y at the level of 
two-point statistics, hence we can express C'' = 4 C" = 4 C^'' (see 
e.g., Bartelmann & Schneider 2001) and subsequently 



(k,x) 



(16) 



In case of the angular galaxy number density power spectrum the 
weight function reads qi = piix) b and subsequently we obtain 



Jo 



b'Ps6{k,x) ■ 



(17) 



We will specify the bias parameter b further in Sect. 5. 



3.2 Projected cross correlation power spectra 

Next we consider the cross correlation power spectra between our 
observables starting with the shear-magnification power spectrum 

v4 rn „2/ 



x)- 



(18) 



JO "^Or) 

The corresponding relation for the shear-galaxy position power 
spectrum reads 

C-(04(^f Q„ r\x'-^brP.Ak,x), (19) 

2 \ c / Jo a(x)x 

where r denotes the correlation parameter for which, similar to the 
bias, we postpone an exact description to Sect. 5. 
Finally, we obtain 

3 Hi r>"' g(y) p(y) 

C™(0=^nn, Ax , brP,,{k,x) (20) 



as the expression for the magnification - galaxy position power 
spectrum. 

We note the following interesting relations, which occur as a con- 
sequence of the polar-scalar filter functions J2 



and 

c''«(^) = ic'He) 



(21) 
(22) 

(23) 
(24) 



We note that these relations can be used to create linear combina- 
tions that can asses the impact of systematics on individual probes. 
We will expand on this in Sect. 5.2. 

3.3 COSEBIs formalism 

The COSEBIs' formalism was developed in SEKIO; we refer the 
reader to this paper for details beyond the brief summary presented 
in this section. 

Throughout this paper we only consider filter functions that are log- 
arithmic in & as these filter functions comprise the second order 
shear information into significantly fewer COSEBI-modes com- 
pared to filter functions that are linear in §. 



3.3.1 Weak lensing 

The COSEBIs shear E-mode, denoted as £„ , can be expressed as 
an integral over the shear 2PCF ^+ as 



2 1„„„ 



(25) 



Note that for a properly constructed TjJ as described below the 
coiTesponding TjZ can be readily calculated as 



XO Id] 7fl2 



(26) 



For further details on this the reader is referred to SEKIO and ref- 
erences therein. In the following we only describe the construction 

of T'^^ 

In order to allow for a proper E/B-modes separation using a 2PCF 
over only a finite interval the filter functions TJJ must meet the 
requirement 



In additic 



(27) 



In addition, the set of filter functions rJJ must be orthonormal, i.e. 



d»Tj,l(§)TZ(i}) = S„, 



(28) 



The explicit construction of the logarithmic T„+ is described in 
SEKIO. The main steps of the construction are: 

• A variable transformation j? — > z = ln(i?/j?mi„) 

• Expressing Eqs. (27, 28) in z with rJJ(!?) tj,l{z) 

• Expanding each tHiz) = 'Z'jto ^'"j^" 

• Calculating the coefficients c„j from the conditions (27, 28) 

Given n, the filter function TJJ will be of order « -1- 1 as it needs 
to fulfill n + I constraints, i.e. it must fulfill Eq. (28) for all TjJ_^_ 
with m < « - 1 and additionally it has to meet the two EB-mode 
separation constraints in Eq. (27). This implies that Tj^ is of order 
two. 



3.3.2 Extension to clustering and magnification 

Having determined the weak lensing COSEBIs filter functions rJJ, 
we can calculate the other five probe's COSEBIs similar to Eq. (25) 



di?j?r,f' 

Jl'min 



(29) 



In this paper we assume the same configuration space filter function 
Til for all other probes, henceforth neglecting the superscripts. 
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Figure 2. Schematic illustration of the modeling of the joint COSEBIs co- 
variance for a given cosmology. 



In addition to Eq. (29) the COSEBIs can be calculated directly from 
the power spectrum 



(30) 



which we only use as a consistency check since computing 
C^^if) ^^^{ff) using a fast Hankel-transformation and subse- 
quently carrying out the finite integration in Eq. (29) is significantly 
faster. 

The Fourier filter functions W^^il) are needed however for the 
computation of the COSEBI's covariance; they can be obtained 
from the T„ as 

lV,7(f) = IVf (f) = IVf (f) = Wf (f) = J" d&&T„(&)Jo(m, (31) 

Wyj{i)=WJ,''if) = J d&&T„(ff)J2m, (32) 
where the yo/2 are a consequence of Eqs. (5 - 10). 



4 MODELING OF COVARIANCES 

We first describe the modeling of covariances for projected power 
spectra; the expression for computing the COSEBIs covariance 
from the power spectrum covariance is straightforward, however 
the actual computation is easily affected by numerical uncertain- 
ties. We outline our method and cross checks at the end of this 
section. 



4.1 Power spectrum covariances 

Under the assumption that the density field is Gaussian (which 
means that the shear four-point function can be expressed in terms 
of two-point functions) the covariance of projected shear power 
spectra can be expressed as (Hu & Jain 2004) 

Covg (c'"'(f)C'^"(f )) = (AC'^''(^)AC™(f )) 
2n 6c,e' 




Figure 3. Full (Gaussian+Non-Gaussian) COSEBIs correlation matrix. 
Since we assume five modes for each of the six probes our data vector 
contains 30 data points, hence the covariance is 30 X 30. We indicate the 
corresponding auto-covariance block matrices in the plots. 



where the superscripts are to be replaced with y, g, fi depending on 
the probe under consideration and £1, denotes the survey volume. 
The covariance gets contributions from the signal C{C) and a noise 
term A'. Note that 



2«=al 



A',,, 



"gal 



1 

"gill 



(34) 



and all other noise terms are zero. We assume the intrinsic shape 
noise, cr^ = 0.32, and note that the factor "2" in the denominator 
results from the fact that the shear has two components. For the 
magnification noise parameter we follow the arguments in Krause 
et al. (2013) defining cr^ = 2cr,//''' for scaling relation based 
estimators of magnification (Huff & Graves 201 1), where / denotes 
the fraction of galaxies for which magnification is measured and cr^ 
being the scatter of the convergence estimator. Being optimistic that 
the method in Huff & Graves (2011) can be extended to late type 
galaxies we assume cr^ = 1.2, noting that the uncertainty of this 
noise level is large. 

Since non-linear structure growth at late time induces significant 
non-Gaussianities in the density field, Eq. (33) underestimates the 
error on cosmological parameters and needs to be amended by an 
additional term, i.e. Cov = Covg + Covng. We model non-Gaussian 
covariance as the sum of the trispectrum contributions (Cooray & 
Hu 2001; Takada & Jain 2009), including a sample variance term 
which describes the scatter in power spectrum measurements due to 
large scale density modes (Takada & Bridle 2007; Sato et al. 2009), 



Covng(C'' 



'(/l),C^°(/2)) 



= -f 

J|l|l 



£1' 
T 

Aih) 



'(1,-1,1,-1) 



(33) 



(35) 
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with r-^'^^a, -1,1', -1') defined as 



Table 1. Cosmological parameter ranges used in tlie likelihood analyses. 



8( X)X 

Or) 



X {p(x)br"-^ x-'T'''4^-^,^^,x\, (36) 
\x X 



where we assume the correlation parameter r = 1 and tt,/? 6 [0; 4]. 
For example the pure shear tri-spectrum T^y^'' , and the pure galaxy 
position tri-spectrum T^sss j-gad 



T'\luh) = 
respectively. 



Jo 



X \x X 



,(37) 
(38) 



4.2 Halo Model Trispectrum 

We model the dark matter trispectrum using the halo model (Sel- 
jak 2000; Cooray & Sheth 2002), which assumes that all matter is 
bound in virialized structures that are modeled as biased tracers of 
the density field. Within this model the statistics of the density field 
can be described by the dark matter distribution within halos on 
small scales, and is dominated by the clustering properties of halos 
and their abundance on large scales. In this model, the trispectrum 
splits into five temis describing the 4-point correlation within one 
halo (the one-halo term r"'), between 2 to 4 halos (two-, three-, 
four-halo term), and a so-called halo sample variance term r^sv^ 
caused by fluctuations in the number of massive halos within the 
survey area. 



T - Til, -I- (r2h,(22) + 7'2h,(13)) + T'sh + 7'4h + T^^^ 



(39) 



The two-halo term is split into two parts, representing correlations 
between two or three points in the first halo and two or one point 
in the second halo. As halos are the building blocks of the density 
field in the halo approach, we need to choose models for their in- 
ternal structure, abundance and clustering in order to build a model 
for the trispectrum. Our implementation of the one-, two- and four- 
halo term contributions to the matter trispectrum follow (Cooray & 
Hu 2001), and we neglect the three-halo term as it is subdominant 
compared to the other terms as the scales of interest for this analy- 
sis. Specifically, we assume NFW halo profiles (Navarro, Frenk & 
White 1997) with the Bullock et al. (2001) fitting formula for the 
halo mass-concentration relation c(M, z), and the Sheth & Tormen 
(1999) fit functions for the halo mass function and linear halo 
bias b(M), neglecting terms involving higher order halo biasing. 



parameter 


3D analyses 


5D analyses 


fiducial 




[0.18;0.38] 


[0.24; 0.29] 


0.2646 


o-g 


[0.65; 0.95] 


[0.76; 0.84] 


0.801 


Wo 


[-0.5; -2.0] 


[-0.8; -1.2] 


-1.0 


ris 


0.963 


[0.85; 1.06] 


0.963 


Wa 


0.0 


[-2.0; 2.0] 


0.0 


b 




[0.9; 1.3] 


1.2 


r 




[0.8; 1.1] 


1.0 



4.3 COSEBIs covariances 

Following Sato et al. (2009) the halo sample variance term can be 
is calculated as 



C0VHSv(C'^''(/l),C™(/2)) = 1^ 



-a 



3^. 

2 £-2 

dxd£l] \ aix) ) ' 



Jo 

X j dM^b(M)l^jj \u(h/x,c(M, 
X j dM'-^fc(M')(yj \u(h/x,c(M' 

X f ^PfiKz{x)W{kx®,)\\ 
Jo 2;r 



zixW 



z(x)f 
(40) 



with u(l[ Ix, c(M, z(x)) the normalized Fourier transform of the halo 
density profile. 

Adding Gaussian (Eq. 33) and non Gaussian covariance (Eq. 35) 
and subsequently integrating over £ and {' we obtain the final 
COSEBI's covariance 

Cov(£f,£™) = ^ df mr(/)iy„'J"(/)CovG(/) 

+ r defwf{£) f df rw™(r)CovNG(/i,/2) 

Jo Jo 



(41) 



This integration must be tested thoroughly for convergence and sta- 
bility of the result with respect to numerical integration precision, 
upper and lower limit of the integration. If the result is stable we 
verify that the covariance and its inverse is positive definite. In 
Fig. 3 we show the correlation matrix of the full covariance ma- 
trix. Since the COSEBIs filter functions average over all Fourier 
modes/angular scales, it is difficult to have an intuitive understand- 
ing of this matrix. We show the impact on likelihood contours when 
neglecting the Non-Gaussian terms in Fig. 5 and when neglecting 
the cross terms of the matrix in Fig. 6. 



5 LIKELIHOOD ANALYSIS 

In this section we simulate various likelihood analyses using the 
data vector in Eq. (1) and covariance in Eq. (4) or subsets thereof 
for survey parameters close to what is expected for the Dark En- 
ergy Survey (flj = 5000 deg^,n„.^[ = 10/arcmin^). The range of the 
cosmological parameter space considered in the different analyses 
is summarized in Table 1. We calculate the likelihood as described 
in Eq. (2) and subsequently the posterior probability via Eq. (3). 
All posteriors are computed on a three- and five-dimensional hy- 
percube, respectively, and then projected onto two dimensions. The 
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Figure 4. Likelihood analysis in three-dimensional cosmological parameter space {left: Sl^ vs trg marginalized over wq, right: Q.^ vs vv'o, marginalized 
over erg). We show the 68% and 95% credible regions for six different data vectors, the solid lines correspond to constraints when using the single probes 
individually, the dashed lines result from a joint analysis of shear and clustering data (red) E = (E'''', E''*, E^*), and when using the full six second-order 
statistics from clustering, shear, and magnification (black) E = (E'''', W, E^*, E^'', E''*, E''^) 



contour plots in Figs. 4 -7 show the marginalized probabihty calcu- 
lated as 

Jr.pJ„(ma.v) 
dPcoi(EIPco) 

where p^^, denotes the remaining cosmological parameters when 
subtracting the considered two parameters from the full parameter 
set. 

The schematic illustrations in Fig. 1 and Fig. 2 show the computa- 
tion of the data vector and covariance, respectively. Note that the 
COSEBIs data vector is computed from the corresponding coiTela- 
tion function over an interval of [r;400']. Computing time for the 
full COSEBIs data vector is < 1 sec per cosmology. 

Comparing individual and combined probes In our first likeli- 
hood analysis we compare the cosmological information of indi- 
vidual probes to a CP analysis. Figure 4 shows constraints from 
the single probe data vectors W', E'''', and E-"* (black, red, and 
green contours, respectively) and their corresponding covariance 
matrices (submatrices of Eq. 4) to the six-probe data vector (dark 
blue contours) described in Eq. (1) and a three-probe data vector 
E = (E'''', E''*, E^*) (light blue contours) which contains only in- 
formation from shear and position. In this analysis we only con- 
sider a three dimensional cosmological parameter space (f2n„ erg, 
and Wo) for the reason that the individual probes error bars cover a 
wide range of cosmologies (which broaden even further for higher 
dimensional parameter spaces), whereas the CP contours require a 
high resolution of grid points. The combination of high resolution 
and large parameter space given our limited computing resources 
and brute-force likelihood approach limits this part of the analysis 
to 3D. 

However, already in the 3D case the results are clear: First, there 
is significant improvement for the CP data vector over the individ- 
ual probes and second, the inclusion of magnification as a probe 
at the assumed noise level does not substantially contribute to the 
information in this parameter space. The latter result also holds for 
the five-dimensional likelihood analysis (see Fig. 6) and can be ex- 



plained by the fact that the information is highly degenerate with 
shear. 

Comparing Gaussian and Non-Gausian covariances Compared 
to Fig. 4, we extend the considered cosmological parameter space 
in Fig. 5 by the spectral index of the power spectrum its- We com- 
pare Gaussian and Non-Gaussian covariances and their impact on 
parameter estimation; for likelihood analyses of individual probes 
such comparisons have been earned out in previous papers (see e.g. 
Takada & Jain 2009; Eifler, Schneider & Hartlap 2009, for cos- 
mic shear), this however is the first time that such a comparison 1) 
is shown for the CP case and 2) includes the Halo Sample Vari- 
ance term in the covariance. By comparing the relative change in 
likelihood contours of Fig. 5 and Fig. 7 we find that both effects, 
i.e. neglecting Non-Gaussianity and assuming perfect knowledge 
of galaxy bias, have comparable impact on the constraints. 

Comparing different CP data vectors Here, we further ex- 
tend our analysis by one cosmological dimension and include a 
simple time-dependent dark energy component w„. For a non- 
tomographic analysis as earned out in this paper constraints on 
Wa are rather poor; overall our marginalized constraints will im- 
prove substantially when extending our formalism to tomogra- 
phy. Figure 6 compares two different data vectors, namely E = 
(E''^E'"',E«s,E^'',E''«,E''S) (black and blue contours) and E = 
(E''^, E-"*, E''^) (red contours); in the first case we further compare 
the full covariance including cross terms amongst probes (black) to 
the case where we ignore these cross terms and only consider the 
auto-coiTelation block matrices in the covariance matrix (blue). 
We find a slight but not very significant improvement of the cos- 
mological information when including the magnification terms in 
the data vector, for the same reasons as outlined in Sect. 5. Using 
only the block diagonal auto-correlation covariance matrix we see 
a notable change in the contours, in particular we also find a change 
in the degeneracies within the cosmological parameter space. Also 
here, we compare the relative change in contour size in Fig. 6 (blue 
and black) to the increase in contour size when dropping the as- 
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Figure 5. Likelihood analysis in four-dimensional cosmological parameter space. We show the 68% and 95% credible regions and marginalize over the all 
other parameters not shown in a given panel. The colors coiTespond to different data vectors and covariances used in the likelihood analysis: black: E = (E'''', 
Wi', W^, W, W^, W^) data vector and full Non-Gaussian covariance, red: E = (E'''', E''''. E^*^, E''''. E''*, W^) data vector but using the Gaussian covariance 
instead. 



sumption of perfectly known bias (see Fig. 7). We find that ne- 
glecting the cross terms has a comparable impact on the change in 
contour size as assuming perfect knowledge of galaxy bias. 



5.1 Uncertainties from Bias and Correlation Parameters 

Understanding the relation of galaxies and their dark matter en- 
vironment is an important aspect of any cosmological parameter 
estimation that includes clustering information. Constraining and 
modeling this relation is an active field of research in theory (e.g., 
Zheng et al. 2005; McDonald & Roy 2009) and observations (e.g., 
Cacciato et al. 2012; Mandelbaum et al. 2012; Julio et al. 2012). In 
practice, any bias model will have to be finetuned to the considered 
data set (galaxy population/morphology and redshift distribution). 
Guidance on any parametrization from first physical principles is 
limited; measurements rely mostly on configuration space quanti- 
ties, i.e. a parametrization in r, z or &. 

Bias and correlation parameter can also be parametrized (e.g. 
Bartelmann & Schneider 2001; Bernstein 2009) as a function of 
k, z in Eqs. (17, 19), more precisely 



(42) 



and 



^lPJKx)PJM^ 



(43) 



with Pgg being the observable galaxy number density power spec- 
trum and Psg being the galaxy-dark matter cross power spectrum. 
In contrast, we parametrize b and r directly as a function of the 
quantity that enters the likelihood analysis, i.e. we define the rela- 
tion between the observable galaxy number density COSEBIs 
and the dark matter density COSEBIs as 



■ X„ Et' 



and similarly 



(44) 



(45) 



For the case of scale-independent bias and coiTelation parameter, 
X„ = and Y„ = br for all modes n. 

Scale (and redshift) dependence of b and r will result in different 
values of X„ (for different modes n) with a correlation induced by 
the functional form of the scale dependence (this also holds for Y„ 
and also for the relation between X„ and Y„). Below, we address 
the issue of how to determine these relations further by proposing 
a measurement scheme to constrain X„, Y„. 
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Figure 6. Likelihood analysis in five dimensional cosmological parameter space. We show the 68% and 95% credible regions and marginalize over the all 
other parameters not shown in a given panel. The colors coixespond to different data vectors and covariances used in the likelihood analysis: black: E = (E'''', 
W^', E^^, W^', W^, W^) data vector and full Non-Gaussian covariance, red: E = (E'''', E''^', E**) data vector and corresponding Non-Gaussian covariance, 
blue: E = (E'''', W^', E**, E'''', E''^', W^) data vector but using only the diagonal block matrices in the covariance. 



In this analysis we bracket X,, e ['^^i„;^max] ™d Y,, e 
[^min'min; ^max''max] with the hmits On b,r in Table 1, which were 
conservatively chosen using guidance from observational results 
from the Red-sequence Cluster Survey (Fig. 9 in Cacciato et al. 
2012). These upper and lower limits are bracketing all X„, Y„, even 
if b and r are functions of k, z- 

Figure 7 shows the results of our simulated likelihood analysis in 
flnij Cg, Wo, n„ Wa parameter space assuming perfect knowledge 
of bias and correlation parameter (red contours), and using the 
parametrization outlined in Eqs. (44, 45). Our bias parametrization 
is conservative in the sense that we do not assume any correlation 
between X„ and Y„ other than that upper and lower limits are calcu- 
lated from upper and lower limits of b and r. 
As expected marginalizing over uncertainties in bias and congela- 
tion parameter significantly weakens cosmological constraints, es- 
pecially for CTg. 

In the following we suggest a measurement framework to make 
progress on the brackets of X„ and Y„ observationally. In the con- 
text of the aperture mass dispersion this method has been suggested 
in van Waerbeke (1998) to detect scale dependence of galaxy bias 
by combining second-order statistics of shear and clustering (see 
Hoekstra et al. 2002; Cacciato et al. 2012, for application to data). 



For the COSEBIs the corresponding relations read 

= fx (46) 

and 

Y„ = fy " . (47) 

The functions fx and /,, depend weakly on cosmology (Hoekstra 
et al. 2002; Schneider, Kochanek & Wambsganss 2006), which 
can be mitigated even further by employing strong priors from 
independent experiments, e.g. Planck. This method allows to mea- 
sure and constrain a "mode-dependent" bias for the COSBEBIs. 
We however note that on cosmological scales galaxy bias has little 
scale but significant redshift dependence and that this method 
should be applied within sufficiently small tomography bins. 
Analogous relations can be derived using E'^,'^ instead of Ej/ and/or 
using Ef/ instead of Ejf . We emphasize that magnification can 
provide additional information to constrain the relation between 
dark and luminous matter; at the very least it provides a valuable 
cross check/nulltest for the above method. 
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Figure 7. Likelihood analysis in five-dimensional cosmological parameter space using the E'''', W,W^, data vector and the con'esponding Non- 

Gaussian covaiiance. Black contours con'espond to a likelihood analysis assuming perfect knowledge of bias and red contours correspond to marginalizing 
over 2 bias parameters h e [0.9; 1.3] and r e [0.8; 1.1]. We show the 68% and 95% credible regions and marginalize over the all other parameters not shown 
in a given panel. 



5.2 NuUtests involving shear and magnification 

Tlie fact that including magnification contributes only little to the 
cosmological constraints is not unexpected given the degeneracy of 
shear and magnification. This degeneracy however allows us to test 
for and to quantify systematics that affect both probes differently. 
As an example we will assume that one of the most important con- 
taminations of cosmic shear, i.e. intrinsic alignment, does not affect 
magnification. We can express the observed shear power spectrum 
as a the sum of the true shear power spectrum and the two intrinsic 
alignment components II (correlation of the intrinsic ellipticity with 
the local density field) and GI (coiTelation of foreground galaxy el- 
lipticity with background shear) (Hirata & Seljak 2004; Bernstein 
2009; Joachimi & Bridle 2010) 

CobsW = a''{t) + C^\l) + C°\l), (48) 
Cl(^) = C'"'(f). (49) 

Using Eq. (22) in terms of the COSEBIs we can rewrite Eqs. (48, 
49) as 

+ Ef = E'll'iobs) - 4 Emobs) , (50) 

thereby constraining intrinsic alignment. We note that contamina- 
tions similar to lA for shear might exist for magnification as well; 



correlations of the local density field with the intrinsic size of the 
galaxies (II''), and/or coixelations of the foreground galaxy size 
with the magnification of a background galaxy (GF) are likely. As 
discussed in Schmidt et al. (2012) magnification estimators which 
are not based on the excess of number densities but on size mea- 
surements have little correlation with their environment (e.g. Cro- 
ton et al. 2005; Maltby et al. 2010) therefore potentially allowing 
the above technique to be successful. 

In any case the above nulltest can be used for other types of con- 
taminations which dominate the shear related quantity but do not 
affect magnification; the best example probably being shear cali- 
bration. 

The number of nulltests is not limited to the magnification and 
shear auto-correlations but additional constraints can be gained 
from three other relations similar to Eqs. (21, 23). As a prerequisite 
we define new COSEBIs filter functions for yg and y/j, denoted as 
r,' such that 



wj,^{n = wjso = wj,''{e), (51) 

^ d&&T„(&)Jo(m = J &&&T'„{-ff)Ji{m, (52) 
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with T„ still being the original filter function TJJ defined in Sect. 
3.3.1. Calculating the new COSEBIs Ejf and Ej,'^ as an integral 
over the coixesponding correlation functions and using the 
we derive the relations 

Eff = 2£•''^ (53) 
= 2Ey, (54) 
£f = 2EJ/, (55) 

which do not hold for the correlation functions (see Eqs. 22, 24). 
These relations or linear combinations thereof can be used to define 
nulltests and subsequently constrain astrophysical uncertainties. 



6 CONCLUSIONS 

In this paper we develop a coherent analysis framework to extract 
constraints on cosmology from all second-order statistics that can 
be derived from a galaxy position (g), shear (7), and magnification 
Qi) catalog. 

For the joint data vector E = (£'''',£''«,£««,£'"',£'''',£''«) we 
develop a fast prediction code that computes projected quantities 
quickly from a density power spectram generated by the Coyote 
Universe Emulator or a modified Halofit implementation. We 
also developed code to predict non-Gaussian covariances for all 
aforementioned projected quantities which is based on the halo 
model. We then generate a CP data vector from our prediction 
code assuming a fiducial cosmology and test this data vector (and 
subsets thereof) in several likelihood analyses, the most extensive 
one covering five cosmological dimensions Q.^, erg, wq, «s, Wa and 
a two parameter galaxy bias model. 

The analysis scheme suggested in this paper differs from previous 
work in several ways: First, we include all second-order cross 
statistics of the observables (g, 7, yu) into the data vector (thereby 
increasing the sources of cosmological information), second, we 
model all cross terms in the covariance matrix, and third, we 
include all higher order correlations of the density field in the co- 
variance matrix. Furthermore we employ the COSEBIs formalism 
to quantify the information content, thereby solving the cosmic 
shear E/B-mode problem and introducing a data compression 
scheme for the other five two-point statistics. 
Not surprisingly, we find substantial improvement in parameter 
constraints when using the CP data vector instead of the individual 
probes, and a sizable increase in the CP likelihood contours when 
fitting for galaxy bias instead of assuming it perfectly known. 
The most interesting result of this paper is the change in likelihood 
contours when including the Non-Gaussianity of the density 
field and all cross terms between the individual probes. We find 
that neglecting either of these features has a comparable impact 
on likelihood contours as the uncertainty of galaxy bias. This 
statement of course only applies to the CP case considered in this 
paper; when only considering galaxy clustering, modeling galaxy 
bias will still be the dominant uncertainty. We however emphasize 
that forecasting exercises for a CP analysis similar to ours should 
incorporate Non-Gaussian covariances including all cross terms. 
The inclusion of magnification slightly improves cosmological 
constraints although at the assumed noise level cr^ = 1.2 this 
improvement is small. We note that magnification in contrast to 
shear and clustering is a relatively recent cosmological probe, 
hence our assumed noise level might be too pessimistic. The 
complete degeneracy of shear and magnification however allows 
for interesting constraints on systematics, e.g. intrinsic alignment. 



shear calibration, etc. We outline several relations that hold in 
the absence of these systematics and suggest extensions of these 
nulltests. 

Scientifically interesting is the extension of our framework to 
other second-order statistics whose distribution follow the same 
likelihood function, e.g. CMB polarization, CMB lensing, CMB 
temperature correlations, and of course also the extension to 
tomography. It is straightforward to apply this analysis scheme 
to a any data set from which can measure projected correlation 
functions. Before extracting meaningful information from data 
however, the framework described here needs extensions. For 
example, implementing a detailed HOD-model approach (van den 
Bosch et al. 2012), adding the parametrization of nuisance param- 
eters, such as baryons, intrinsic alignment, photo-z calibration, 
and shear calibration is required. Although this is challenging 
infrastructure work, the inclusion of these add-ons poses no 
fundamental problem for our framework. 
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